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Abstract 

q!^ ■ The dipole coupling term between a system of N particles with total charge zero and the 

' electromagnetic field is derived in the presence of a weak gravitational field. It is shown that 

the form of the coupling remains the same as in flat space-time if it is written with respect 
to the proper time of the observer and to the measurable field components. Some remarks 
O ' concerning the connection between the minimal and the dipole coupling are given. 

o" 
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1 Introduction 



Although modern theories of quantized matter in curved space make predictions for extreme situations like, 
Q\ ' e.g., the very early universe still very few experiments are done to shed light on the connection between 

general relativity and quantum mechanics. The recent progress in atomic interferomctry may lead to new 
contributions on this topic. For example, the Sagnac phase Q and the influence of the earth's acceleration 
were measured showing the behaviour of atoms in non inertial frames of reference. These experiments, 
' which use the interaction between lasers and atoms to split and recombine the atomic beam, may also be 

, of use to measure the influence of space-time curvature on atoms [|| . It is therefore of interest to study the 

""^ ' behaviour of such devices in weak gravitational fields like that of the earth. In order to do so it is necessary 

q-i to generalize the theoretical methods which are used for the description of atoms and lasers from a flat 

jjj^ 1 space-time to a curved one. 

{jjQ, In quantum optics it is often convenient to use the dipole coupling —d-E instead of the minimal coupling 

scheme in the calculations. While the latter is invariant under Lorentz transformations the dipole coupling 
has the advantage that it is directly related to the physical electric field E and not to the gauge dependent 
vector potential A^. The (non relativistic) equivalence of the two approaches was first demonstrated by M. 
Goppert-Mayer [Q . She showed that the classical Lagrangians of the two theories are related by a canonical 
transformation. Power and Zienau |^| have extended her work by use of a unitary transformation in quantum 
theory and derived, in an approximation, a multipolar Hamiltonian where the dipole coupling is only the 
first order interaction term of a multipole series. This transformation was made exact by Woolley (see, e.g., 
Ref. p). All these derivations are formulated in a fixed non-covariant gauge, mostly the Coulomb gauge. 
The generalization to an arbitrary gauge fixing was done by Power and Thirunamachandran |?J by adding a 
total time derivative to the classical Lagrangian and by Woolley ||. 

Despite the formal equivalence of both couplings there has been a long discussion in the literature which 
one is better suited to describe the interaction between matter and light. Lamb || stated that the dipole 
interaction is preferable because the results of certain calculations fit better to the experiment. The minimal 
coupling scheme on the other hand is covariant and connected with a gauge symmetry, facts which seem 
to imply that this coupling is more fundamental. Ackerhalt and Milonni have pointed out that this 
controversy may arise from the necessity to transform also the states if one applies a unitary transformation. 
This leads to the question for which coupling scheme the textbook wavefunctions are the right choice. We 
will argue in this paper that the textbook wavefunctions seem to belong to the dipole coupling. 

The main subject however is to study the modification of the dipole coupling in weak gravitational fields. 
Our strategy is to begin with a general relativistic Lagrangian for a multiparticle system in an electromagnetic 
field, to consider the limit of small velocities, and to perform the Power-Zienau transformation along the 
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lines of Ref. [|7|. This approach has two advantages: one can avoid to use the Coulomb gauge which is 
not covariant, and one circumvents the inclusion of gauge constraints like, e.g., the Gupta-Blculer condition 
which are not common in quantum optics. One disadvantage is that the spin has to be treated separately. 
To the knowledge of the author the dipole coupling in gravitational fields was addressed only by Borde et 
al. ]IT| , |l2| by replacing the space-indices of the non-covariant expression by tetrad indices. In this paper it 
is the aim to derive its structure from an underlying multi-particle theory. 

We use the conventions of Ref. for general relativity, i.e. sgn g^ v — +2. Greek indices run from 
to 3 and latin ones from 1 to 3. Summation is understood whenever an index appears twice. Tetrad indices 
are underlined. We use natural units (h = c = 1) and Heaviside-Lorentz conventions for the electromagnetic 
field. 



2 The derivation of the Hamiltonian 

We begin with the covariant action S of a system of particles and the electromagnetic field, 

S = { C^/^d i x = [ Ldx° 



(1) 



where 



N 

E 



-1/2 



q^A^W - - / d'yy/^gJyjF^F^iy) . (2) 



a labels the particles which travel on the trajectories X? a y The dot denotes the derivative with respect 



to the coordinate time x° = x a 



y°. Following Rcfs. 



§§ 



we perform the limit of small velocities, 



<C 1, and consider only weak gravitational fields, g^ = r/^ + with <C 1. This leads to the 



new Lagrangian 



L = 



2^ ) ~^~ hw + m («) hoi x («) + ~y~ 

+q (a) A {x {a) ) + q (a) A i (x (ta ))x\ a - ) \ -t 



(1 + -h 00 )S k i + hu 



d 3 y< 



(A kfi - A 0tk )(Ai t0 - Aqj) - 



(1 - -^h X \)5 k i + h k i 



(1 + -jhx\)8 k i - h u 



BkBi — 2(^4^0 — Ao^)eijkBjh 



(3) 



ok 



Each factor of in the sum over a has to be taken at the point X( a y Throughout the paper all expressions 
are calculated to first order in only. In the derivation of (|^) we have subtracted the total rest energy 
M = m( Q ) and have defined Bi := eijkAkj where Eijk is the total antisymmetric symbol with £123 = 1. A 
comma denotes the derivative with respect to the following coordinate. Here it is necessary to give a remark 
concerning the position of the indices. In general relativity an index which appears twice has to appear as one 
upper and one lower index. Any expression which does not fulfill this requirement cannot be invariant under 
coordinate transformations. Although the starting point of this calculation was a covariant expression we will 
use extensively differential geometric methods in fiat three-dimensional space. It is therefore convenient to 
switch to a three-space notation where the indices of any three- vector V = {Vi} are lower case indices except 
for all coordinates X( Q ),?/, . . . where, e.g., y = {y 1 }. This can be done without making errors as long as no 
index is moved with the space-time metric, and as long as we do not transform to another coordinate system. 
The resulting Hamiltonian then describes the time evolution of the particles and fields in this coordinate 
system. 

Eq. (||) is the Lagrangian of non relativistic particles moving in weak gravitational field. If one defines 
a Hamiltonian H = px — L one ends up with the minimal coupling scheme. As will be shown below this 
non-covariant way to define H is only justified in the case of time independent gravitational fields. In order 
to derive the dipole approximation of the coupling we follow Ref. R] and add a total time derivative to the 
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Lagrangian: 



with 



L'=L 



dx° 



d 3 yP ■ A 



P{lj) = y2 Q(a)(X( a ) - Rem) \ 8(y - R cm - A(.f( Q ) - R cm )) dX . 

1 N 



(4) 
(5) 
(G) 



is the center of mass of the particles. This is the point where covariance is explicitly broken by adding a 
non-covariant term to L. This step is necessary as the aim, the dipole coupling, is also non-covariant. Terms 
like y/—g or ^fgr, (gs is the determinant of the metric on the hypersurface x° = const.) are not included in 
the integrand of Eq. because this would destroy the derivation if the definition of P is not appropriately 
modified. Such alternative derivations should lead to the same result if measurable quantities are considered. 

It is not difficult to show that the (flat space) polarization field P is related in the sense of distributions 
to the (flat space) magnetization 



M{y) = ^2 q( a ) (x( a ) - R cm ) x (x (q ) - R cm ) / X8(y - R cm - X(x (a ) - R cm )) dX 



and the Rontgen current Jro = rot(P x R cm ) via 

P + mtM + J R6 = ^2 s (y- x (a) ) - QR cm 8(y- R c 



(7) 



(8) 



In the remainder we assume the system to be neutral, Q = J2 a 1{a) = 0- ^ n addition, we decompose vectors 
V — + V 1 - related to the electromagnetic field into its longitudinal (rot V" — 0) and transverse ( div 
V 1 - = 0) part, for example 

A = A 11 + A 1 , rotA 11 = divA 1 = . (9) 

Note that the spatial integral over the scalar product of any transverse vector with any longitudinal vector 
vanishes and that the longitudinal part of the polarization for Q — is given by 



P 11 = Wn 



where 

With Eq. (||) we arrive at 

* = £ 



V Q (x) = -- d 3 y 



-in 



divP(y) _ 1 



-Y 



\x-y\ An ^ \x- x (a) \ 



(10) 
(11) 



— y~ h Q + m (a) h Q ■ x [a) + 

d 3 y< 



(1 + 2^oo)4i + h kt 



(1 + -h\x)8ki -hid 



{^--jh\x)8u + hki 



(Ak,o-A , k )(Ai, -Ao,i)- 
2\A - V. In i • (B < hn) } -!- I d 3 y{l ± ■ [rotA/ + J R5 ] - A^P 1 - - AP 11 j 



B k Bi 



(12) 



For notational convenience we have introduced the vector /i = {^oi}- The new Lagrangian L 1 has the feature 
that it depends only on Aj , not on A| which is therefore a cyclic variable. The corresponding Routhian (see, 
e.g., Ref. m&) is given by 



R = L' 



d 3 yA 11 



dA 



ill 



(13) 
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The derivative in the r.h.s. is given by 



which is a constant of motion. The vector A is defined by 

1 



A; 



1 + 2^ AA ) _ 



(A 3 ,o - A j) + (P x /lo^ + Pn 



(14) 



(15) 



and agrees in absence of a gravitational field with the electric displacement. Using Eq. (|T^) to eliminate the 
cyclic variables from R one finds an expression with no time derivative of Aq and which depends linearly 
on it. Hence, Aq plays the role of a Lagrangian multiplier. Solving the corresponding constraint leads to 
pi — A.!! = VVb which implies A" = 0. We thus find for the Routhian 



R = 



E 



m(a) r ■-. , m (a) 

——/lOO + "1(a) flo • X( a ) + — — 



(1 + -h 00 )S k i + hki 



*(a)*(a) f _ ^o«i 



(16) 



<i 3 ?/< (A ) 2 - 2P^ ■ ^4 + 2P-(M + PxP) 



(1- -jh X \)5ki + hki 



BkBi 



-2(PH + A ) • (P x ho) + 



(pn+4 ) fc (p"+i ), 



where 



8?r ^ 1^ 



9a 9/3 



= 1 / d»y(J>ll)> 



(17) 

Qlj 8 " X '(/3)l ' ^ 

is the total Coulomb interaction between all charges. One may interpret this explicit occurence of the 
Coulomb interaction between the particles as an indication that the textbook wavefunctions for atomic 
electrons, which are usually derived by assuming the interaction to be of the Coulomb-type, belong to the 
dipole coupling. 

The Hamiltonian is defined by the usual relation 



H = J2 



P(a) ■ E(a) 



a = l 



d 3 yn x -A -R 



(18) 



where p/ a \ and H 1 - are the canonical momenta of the particles and the transverse electromagnetic field, 
respectively. Performing the calculations and neglecting terms of the order hfi V /m a we arrive at 



JV 



H = 



a=l 
1 



E TO « ( 1 _ \ h oa(x( a ))) + E 2~~ ( P( a ) ~ m aho{x( a )) - I d 3 yB 



a=l 



d 3 y< 



{^--^hxx)S k i + h k i 



(A^-A/- + B k Bi) - 2A k L P l + P k P t 



-2(P- A x ) ■ (B x h 



(19) 



where the vectors ri a are defined by 



n a {y) 



M 



Jo 



(0) 



idX 



+ q a r( a ) / XS(y- R cm - Xf {a) )dX . (20) 
Jo 



In the last equation we have introduced the nonrelativistic relative coordinates and the center of mass 
position 



(a) ■= x (a) 



(21) 
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which will be of use in the derivation of the dipole approximation in the following section. 

Note that this Hamiltonian contains explicitly the (modified) Coulomb potential. This can be seen 
by using Eq. (|l7]) for the evaluation of the term proportional to Pj^Pi in Eq. ([l9|); we will discuss the 
modifications to the Coulomb potential below. The term proportional to A^P; describes the modified 
dipole coupling between the transverse electric displacement and the total dipole momentum of the atom. 



3 The dipole approximation 

To gain more physical insight into the result (|l^) it is of advantage to perform the dipole approximation 
which is valid for many quantum optical applications. In our context it amounts in the restriction to terms 
which are at most linear in the relative coordinates fv a ) . In this limit the polarization and the magnetization 
become 

P(y) = d5(y~ R cm ) , M(y) = x ^<Z(a) ^O) x ^(«) 5 (v ~ ^m) (22) 



A' 



where 

d := ^ q a r( a ) (23) 

a=l 

is the dipole operator of the atom. It is convenient to exploit the fact that the Hamiltonian transforms as 
a scalar under a change of the dynamical variables. We take as new variables of the particles the center of 
mass position R cm and the relative coordinates f( a ) of the first N — 1 particles. Denoting the corresponding 
canonical momenta as P cm and fir a \ it is not difficult to show that these are related to the old momenta by 



= "T"F I cm ~ P(/3) I P(a) > a = l,...,N- 




Inserting this into Eq. ( |l9|) and performing the dipole approximation we can write the Hamiltonian as the 
sum 

H = H at + Hem + Hro + iJ ra d + -ffint (25) 

of five parts. The first is the internal Hamiltonian 

N-l ^ m N-1 2 N-1 



Hat — 



2m (a) \^ a > M 

a=l /3=1 



+ \l<f i !>{lP-)*- 



\h xx S kl -h kl 



PkPi-+PkPi+2PkPn\ (26) 



which does not depend on the transversal electromagnetic field and P cm . The first sum is the kinetic energy 
term in which the subtraction of the sum over (3 describes the generalization of the reduced mass for more 
than two particles. The last term in the second line contains modifications to the Coulomb potential and to 
the dipole energy. It should be noted that the Coulomb potential contains the self-energy of the particles 
and is therefore divergent. This is also the case for the term proportional to (P^) 2 which describes the 
dipole self-energy. The last sum in the first line describes the coupling of the internal angular momentum 
to a rotation. This can be seen by switching to the Fermi coordinates of an rotating observer in the weakly 
curved space p7| . In this coordinate system the components hoi(x) are essentially given by EijkUJ-'x k where 



is the angular velocity of the observer. Inserting this into Eq. (26) shows that the corresponding term is 



of the form —uj ■ J2 a L( a ) if £■(» is the orbital angular momentum of particle a. 
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The center of mass contribution 



—P 2 

2M cm 



M 



hoo{R c 



Pcm • ho(Rcm) 



(27) 



has the same structure as the Hamiltonian for a free particle in a weakly curved space, comp. Refs. |14[ [if 
The internal and external degrees of freedom of the atom are coupled via the Rontgen term 



H Ro = -^Pcm ■ (B X d) . 



(28) 



This is the same expression as in Minkowski space because we have neglected all terms of the order h^ v /m^ a y 
The radiative part of H is found to be 



ffrad = ^ / <*V 



1 - ~h\\ j <">/,/ + 1 1 1:1 



(A£Aj- + B k Bi) - 2h ■ (A 1 - x B) 



(29) 



and contains a coupling between the Poynting vector and the rotation. Here we have used the relation 

fL x = —K x . 

The interaction between matter and radiation is described by 



1 



h\\ 8m + h k 



A£di - B 



y ■ 



q(a) 



P(a) 



1(a) 



N-l 



AI 



0=1 




(30) 



The first term describes the dipole coupling between the atom and the transverse electromagnetic field. 
The second term is the well known coupling between the magnetic field and the angular momentum of the 
particles. As in Eq. (26) the mass reduction has to be taken into account. Again this term is the same as in 
flat space since we have neglected terms of the order of Oih^y /to( q )). 



4 The dipole coupling between measurable quantities 

To make contact with the measurable quantities of the theory it is necessary in general relativity to consider 
the components of each tensor with respect to a tetrad field which fulfills e^ep^ = r\ap at each point in 
space. Here T]afi is the Minkowski metric. The measured components of, e.g., the electric field are then given 

by 

E l -=F^= ePft&F^ (31) 

(Tetrad indices are raised and lowered with the Minkowski metric). We now focus on the dipole term in 
Eq. (|30|). The vector eo of the tetrad is assumed to be orthogonal to the hypersurfaces x° =const. and is 
therefore given by e[] = 1 + /ioo/2, the rest of its components vanish. The orthogonality of the tetrad then 
implies e;o = for all three space- like vectors e^. 

Consider now the dipole term in Eq. (30) and insert Eq. (p"5|). The dipole coupling then has the form 



-Fqi + {Bx ho)i + 



1 - j iS u 



h 



1:1 



Pi- 



rn) 



Since the longitudinal part of the vector A is zero we can omit in this expression the index _L. First we 
examine the term proportional to Fqi . Recalling that the dipole moment is the sum of relative coordinates^] 



(«) 



it is not difficult to find 



F i % 



'(a) 



(1 



hoo/2)Fm2_^q {a) r- ia) 



= -(1 - h m /2)Ei^q (a) r- {a) 



(33) 



3 To achieve true invariant quantities it would be necessary to work with the derivative d ll (s 2 ^) /2 of the geodesic distance 
S( a ) between the center of mass of the atom and the particle a instead of the relative coordinate of the particle. But 

> 9""(ficm) d„(s 2 (a) )/2 so that the 



since we are working in the dipole approximation (up to linear terms in r^ a )) we have 
difference is unimportant. 



G 



If we interpret 



PF := (1 + h m /2) 



1 - 2 hxx i 



P k + (Bx h ) 



(34) 



as the actual polarisation field modified by the gravitational held we hnd for the dipole term the expression 



H dip = -(l-h 00 /2)D^d k 



(35) 



Compared to flat space the dipole coupling seems to be modified by the factor (1 — hoo/2) which is identical 
to \/—goo- But even this term disappears if we take into account that the Hamiltonian (|2{]) describes the 
time evolution of the system with respect to the coordinate time x° . If we switch to the proper time r of a 
family of observers with four-velocity — the Hamiltonian density H., defined by H = J H.d 3 y, has to 
be transformed according to 

H' = H^r- = Kel = H^L= . (36) 
dr - v-goo 

This implies that the dipole coupling between the atoms and the electric field takes the form 

H' dip = -DiS- (37) 

even in a weakly curved space when it is expressed with respect to the measured quantities and the proper 
time of an observer. This result is in agreement with a recent work of Lammerzahl fl9|| in which he derives 
the dipole coupling for a quantum mechanical particle moving in a PPN space-time. 

We now discuss briefly the modifications to the Coulomb potential as this is of interest for the decision 
to which form of coupling the textbook wavefunctions belong (see the remark below Eq. (|l7j)). It is of 
advantage to start with the expression 



0--^h\x)S k i + h k i 



P k Pi 



(38) 



in Eq. (|19|) which contains the Coulomb as well as the dipole energy of the atom. In a first step we replace 
the flat space polarization P[ by its gravitational counterpart PF and arive at 



V = 



\ J d 3 y^(l-h 



(l+2^AA)4i - h kl 



2P G ■ (B x h ) 



(39) 



The last term can be interpreted as a gravitational induced interaction between the atom and the magnetic 
field and is therefore not a part of the Coulomb or dipole energy; we will omit this part. Performing the 
same steps as for the dipole coupling we find that the first term has, after the insertion of the tetrad vectors, 
the form 

' f d 3 y(l + h x x /2)PFP Gl - . (40) 

This 



It only left to replace the coordinate time by the proper time of the observers according to Eq. 
leads to the result that the sum of Coulomb and dipole energy is given by 



dVPFP 



G u Gi 



(41) 



where dV = (1 + hu/2)d 3 y = ^/g^d 3 y is the three- volume element of the hypersurface a: = constant. It has 
therefore the same form as in flat space if it is written with respect to the measurable quantities and with 
respect to the proper time. 

We thus find the result that in connection with the dipole coupling the modifications of the Coulomb 
potential are related to the modifications of the (longitudinal part of) the polarization field. If the calculation 
of the modifications in the minimal coupling scheme are different from the present result this would, in 
principle, open a new way to test experimentally to which coupling the textbook wavefunctions belong. It 
should be noted, however, that the magnitude of the modifications in atomic systems is far too small to be 
measurable. 
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This section will be closed with some general remarks concerning the use of the family of observers to 
define the tetrads and therefore the measurable quantities. In principle it would be more convenient to 
work in the atom's frame of reference instead of including the family of obervers in the description. But 
in general relativistic situations a frame of reference cannot be unambiguously defined. Usually one uses 
Fermi coordinates pC[ | as a local approximation for the reference frame, but this construction has some 
shortcomings. A modification of it may circumvent these problems |2l]], but the true form of a reference 
frame, if it exists, has to be determined by experiments. 

A second remark concerns the fact that the Hamiltonian is changed if we switch to the proper time of 
the family of observers. This is due to the fact that the Hamilton operator is the time evolution operator 
of the system. If we change the time coordinate then the evolution operator is also changed. The use of a 
proper time is of advantage compared to the coordinate time since it is directly measurable by the clocks of 
the observers. 



5 The spin interaction 



In the derivation of the Hamiltonian the inclusion of the spin was not possible because only classical particles 
were considered. But the heart of the Power-Zienau transformation, the subtraction of a total time derivative 
from the Lagrangian, can be made without any reference to the spin. Hence, it should be possible to include 
the spin by deriving the minimal coupling Hamiltonian and by taking over the spin terms into the Hamiltonian 
with dipole coupling. In order to get the correct spin terms we will follow closely the approach of Ref . || . 
For brevity we will only scetch the main steps and refer to this paper for further details. 
The Dirac equation in a weak gravitational field can be written as 



idotp = Hip 



(42) 



with 



H 



— hQjSjik^k ' 



XV l,p 



h", 



2 ^ij a i 
~2 h °° 



(-idi - qAi) 



7o - .^oiTi 



(43) 



In order to give the scalar product between spinors the usual form in flat space we redefine the field by 



4, 



1 



1 



1 



(44) 



Oib' with O 

y 4 4 

The corrections to the Pauli equation can be found by performing a Foldy-Wouthuysen transformation (see 
e.g., Ref. JL8[ ) with the unitary operator exp(zS') where the Hermitean operator 5" is given by 



S = — 7i(-«9 l - qAi) - - — h ik ^k(-idi 
2m ~ 4m 

I ■ - Hamiltonian for the Schrodinger field is then found to be 



qAi) 



~lihki,k 



Zm 



idi - qAi - mh 0i ) 2 - qA Q + 



1 - 2^00 



-^hoi^sukak 



-Z—D k Ok 

Zm 



This is the same result as in Refs. |TJ, except for the inclusion of the spin-gravity interaction 

— -rhoiiSiikcrk ■ 



(45) 



(46) 



(47) 



Here Ok arc the Pauli matrices of the particle under consideration, say particle a. In the discussion of the 
result ( p5| ) we have seen that the term h i^ is essentially given by ei mi Lo m with the angular velocity uj m of 
the observer. It is therefore obvious that this term describes the — Co ■ S coupling of the spin part S of the 
total angular momentum to the rotation Co. 

The derivation along the lines of Ref. ||] shows also that the Hamiltonian is only related to a relativistic 
Hermitean energy operator if the gravitational field is time independent. Since the spinless part of it is 
exactly the same as the canonical Hamiltonian, and since the canonical procedure is not manifest covariant, 
we conclude that the same condition must hold for the Hamiltonian (|2£ 



8 



6 Discussion 



In this paper we have shown that the interaction between atomic particles and the electromagnetic field 
can be described with a dipole coupling term also if space-time is weakly curved. The form of the coupling 
remains the same as in Minkowski space if it is written with respect to the proper time of an observer and 
the measurable quantities of the theory. The same is true for the Coulomb and the dipole energy of the 
atom. 

The central assumptions in this derivation are the smallness of the velocity of each particle, the weakness 
of the gravitational field, and the validity of the dipole approximation. A further assumption was tacitly 
made by using the non-relativistic form of the center-of-mass coordinates. Although this should in general 
be a good approximation, Fischbach et. al. |2^] have shown that the use of different relative coordinates 
(center-of-energy, e.g.) can result in different perturbational contributions of the weak gravitational field 
which do not vanish even when the mass of the nucleus is very large. We assume that these contributions 
are small enough to be neglected. 

It is interesting to make a comparison of the present results with the well known formal equivalence 
between the Maxwell- field in a curved space and a dielectric medium (^3). In this approach one defines a 
formal dielectric displacement vector to describe the influence of gravity on the Maxwell field. In absence of 
particles, i.e. for vanishing polarisation P, the formal electric displacement agrees with the vector A defined 
above (in Ref. [2^] the presence of charged matter was not considered). Also the coupling of the Poynting 
vector to the rotation occurs in the energy density of the formal Maxwell field. 
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